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ABSTRACT 

Context. The exact nature of the quiet solar chromosphere, and especially its temporal variations, are still subject of intense debate. 
One of the contentious issues is the possible role of magnetic field in structuring the quieter solar regions. 

Aims. We characterize the dynamics of the quiet inter-network chromosphere by studying the occurrence of acoustic shocks and their 
relation with the concomitant photospheric structure and dynamics, including small scale magnetic structures. 
Methods. We analyze a comprehensive data set that includes high resolution chromospheric (Call 854.2 nm) and photospheric 
(Fel 709.0 nm) spectra obtained with the IBIS imaging spectrometer in two quiet-Sun regions. This is complemented by high- 
resolution sequences of MDI magnetograms of the same targets. From the chromospheric spectra we identify the spatio- temporal 
occurrence of the acoustic shocks. We compare it with the photospheric dynamics by means of both Fourier and wavelet analysis, and 
study the influence of magnetic structures on the phenomenon. 

Results. Mid-chromospheric shocks occur within the general chromospheric dynamics pattern of acoustic waves propagating from 
the photosphere. In particular, they appear as a response to underlying powerful photospheric motions at periodicities nearing the 
acoustic cut-off, consistent with 1-D hydrodynamical modeling. However, their spatial distribution within the supergranular cells is 
highly dependent on the local magnetic topology, both at the network and internetwork scale. We find that large portions of the inter- 
network regions undergo very few shocks, as "shadowed" by the horizontal component of the magnetic field. The latter is betrayed by 
the presence of chromospheric fibrils, observed in the core of the Call line as slanted structures with distinct dynamical properties. 
The shadow mechanism appears to operate also on the very small scales of inter-network magnetic elements, and provides for a very 
pervasive influence of the magnetic field even in the quietest region analyzed. 

Conclusions. The magnetic field might play a larger role in structuring the quiet solar chromosphere than normally assumed. The 
presence of fibrils highlights a clear disconnection between the photospheric dynamics and the response of the geometrically overlay- 
ing chromosphere. As these results hold for a mid-chromospheric indicator such as the Call 854.2 line, it is expected that diagnostics 
formed in higher layers, such as UV lines and continua, will be affected to a larger extent by the presence of magnetic fields, even in 
quiet regions. This is of relevance for the chromospheric models that make use of such diagnostics. 

Key words. Sun: chromosphere — Sun: magnetic fields — Sun: oscillations 
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, 1. Introduction A landmark work on thi s topic was the 1-D radiative- 

hydr odynamical modeling of Carls son & Stein! ( 1 19951 1 1 9971 

Observations have long shown that the outer solar atmosphere |20Q2|) , who derived the chromospheric response to acoustic 

is not in radiative equilibrium, with the net radiative loss of waves in the absence of magnetic fields. Using an observed pho- 

the quiet chromosphe re estimated at bet ween « 4 kWm" 2 tospheric piston, they could reproduce the temporal evolution of 

(VernazzaetalJ |1981[) and 14 kWirT 2 (|Anderson & Athayl the chromospheric Ca II H spectral line over long intervals with 

119891) . The identification of the source and exact mechanism of remarkable accuracy. In particular, they reproduced the occur- 

deposition of the energy necessary to maintain a stationary situa- renC e of the U 2 v (K 2V ) grains observed in the spectral profile of 

tion is often referred to as the problem of chromospheric heating, the resonance H and K lines of of Ca II. The grains consist of 

It has been a vexing problem in solar physics for decades. quasi-periodic brightenings in the violet wings of the lines, with 

In stellar atmospheres, strong density stratification makes a -150-200 s repetition rate, and ar e associated with strong red- 



it relatively easy for propagating acoustic waves, excited by shifts of the line core position ( see iRutten & Uitenbroek 1991, 



the turbulent convection, to develop into shocks. Indeed this hereafter RU91, and Section IBl. lCarlsson & Steinl dl997l). here- 

mechanism was proposed early on as a plausible means to inafter CS97, showed how the grains are due to weak acous- 

provide the nece s sary energy input to the solar chromosphere tic shocks developing at heights of about 1 Mm in the non 

dBiermannl 1 1 948b ISchwarzschildl 1 1 9481) . Much theoretical and magnetic atmosphere, driven primarily by photospheric acoustic 

numerical work has since been conducted to assess the vi- waves with periodicities close to the cutoff value of -5.5 mHz. 

ability of the acoustic heating as the basic heating mecha- A controversial by-product of their model is that the intermittent 

nism for nonmagnetic chromosphe r es of slowly rotating stars acoustic shocks greatly enhance both the emissivity and the tem- 

(e.g. jNarain & Ulmschneidetj 119961: lUlmschneider & Musielakl perature but only locally, without causing a net outward increase 



l2Q03h ! in the average chromospheric temperature dCarlsson & Steinl 
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1 1995b . If the CS97 shocks correctly represent the acoustic heat- 
ing mechanism, the non-magnetic chromosphere would be es- 
sentially cold. This is contrary to the common view, derived 
from semi-empirical static models, of a chrom osphere charac- 



terized everywhere by a steady temperature rise ( Vernazza et all 
1 198 It iFontenla et all 120071; lAvrett & LoesenbOOaT The failure 



of the CS97 model to reproduce spatially and temporally re- 
solved quiet Sun U V observations ob tained with SUMER in 
chromospheric lines (Carlsso n et al.ll997h has further sharpened 
the controversy. 

A strong debate has subsequently ensued about the pos- 
sible shortcomings of the CS97 model and the general vi- 
ability of the acoustic heating mechanism in shaping the 
quiet solar chromosphere. Much discussion has been de- 
voted to the role of high frequency acoustic waves (among 
many others. | Fossum & Carlssonl2 005: Wed emever-Bohm et al.1 
2007; ICuntz et al.l 120071; ICarls son et al.1 120071) the n eed to 
address the problem in 3-D (e.g. Wedeme yer et al.l [2004; 
lUlmschneider et aTl 120051) , and the r o le of magnetic fields 
in th e quiet S u n (e. g. Pudge et all 120031: iMcIntosh et alj 
120031: ICarlssonl 120071; IWedemever-Bohm et al.1 120071) . At 
the bottom of the debate lies the fundamental question 
of whether the internetwork chromosphere is wholly dy- 
namic in natur e (Carlsson 2007; Martfnez -Svkora et al.l 12008; 
IWedemeyer-Bohm & Wogerl 2008) or whether the dynamic vari- 
ations represent only minor perturbations on a s emi-static state 
of the kind described in sem i-empirical models (Kal kofen et al.l 
1 1999t lAvrett & Loeseril2Tj08l 

From the observational point of view, gaining further insight 
into this issue requires new data that can simultaneously achieve 
high spatial and temporal resolution, and spectral information 
with different diagnostics, over extended fields of view (FOV). 
Most of the debate to this point has been based on (older) spec- 
trographic observations obtained in the Ca II H and K lines, and 
on UV spectral data obtai ned by Skylab or, m ore recently, by 
SUMER on board SOHO dWilhelm et al.lll995l) . Spectrographic 
data suffer however from a limited FOV in fixed slit observa- 
tions, or from a low cadence when scanning an extended region. 
In the UV, the limited flux and instrumental response further 
lower the achievable spatial and temporal resolution, an impor- 
tant limitation if the dynamic picture just described for the chro - 
mosphere is the correct one (Wedemeve r-Bohm & Wogerl2 008). 
Historically, imaging with broad-band filters has been widely 
used, both in the UV and visible range. Such data, however, of- 
ten mix signals arising from vastly different regions of the so- 
lar atmosphere, as well as miss signatures of strong velocities 
that might shift the lines outside the filter band. For a general 
review on the issues of chromospheric observations we refer to 
lRuttenl(|2007l). only add ing here that recently the POLIS spectro- 
graph (Be ck et aT1 l2005) has become operational, producing new 
high quality Ca II H spectr a that have been employed in studies 
of the quiet chrom osphere dRezaei et alJl2007llBeck et al.f 2008; 
iRezaei et al.ll2008l) . 

In this paper we use an entirely different type of in strument, 
the In terferometric BIdimensional Spectrometer (IBIS, Cavallini 
2006), to derive novel results about the dynamics of the quiet 
solar chromosphere and in particular the presence (or absence) 
of acoustic shocks. IBIS is an imaging spectrometer installed 
at the Dunn Solar Telescope (DST) of the US National Solar 
Observatory, and provides observations that combine high spa- 
tial and temporal resolution over an extended two-dimensional 
FOV, with the full spectral information in both photospheric 
and chromospheric lines. In this sense, it is an ideal instrument 
to overcome many of the observational shortcomings outlined 



Fig. 2. Full disk EIT images acquired in Fe XII 19.5 nm on 2004 
May 31 (left) and 2004 June 02 (right). The cross-hairs indicate 
the center of the IBIS FOV for each day. The IBIS FOV covers 
roughly the small coronal bright point intersected by the cross- 
hair on the right panel. 



above. We analyze the behavi or of the chromosp heric Ca II 854.2 
nm line which, as shown in Cauzz fet alj (120081) (hereafter Paper 
I), is one of the most promising diagnostics for high resolu- 
tion chromospheric studies. Equally important, the availability 
of spectral information over a 2-D FOV is crucial to understand- 
ing the very dynamic and non-local chromosphere. We further 
combine the Call 854.2 observations with simultaneous and 
cospatial photospheric data also obtained w ith IBIS and mag- 
netic data from MDI dScherrer et alj[l995l) . This allows us to 
address the relationships between the chromospheric and photo- 
spheric dynamics, and the influence of magnetic field. We an- 
alyze two completely analogous data sets obtained on separate 
days, both times in quiet regions near disk center, and in so do- 
ing uncover the fundamental role played by the local magnetic 
topology. 

The paper is organized as follows. In Section [2] we de- 
scribe the observations. Sections [3] and [4] provide evidence that 
mid-chromospheric acoustic shocks are clearly observed in the 
Ca II 854.2 nm line, and outline the methods we use to identify 
them. In Sect. we report on the derived shock properties, both 
in relation to earlier findings obtained primarily with Call K 
observations, and shockl5.gc.as new results regarding their spa- 
tial distribution. Sections [6] and [7] address the occurrence of the 
shocks in relation to the magnetic field structure and evolution, 
and the photospheric dynamics. Finally, in Sections [8] and [9] we 
discuss our findings and provide conclusions. 



2. Observations 

The IBIS design and gener al issues on data reduction have been 
described in earlier papers (ICavallinil2006l ; lReardon & Cavallinil 
120081; |J anssen & Cauzzil 120061) . so we mention here only the 
characteristics most relevant to the present work. 

IBIS is based on two air spaced Fabry Perot Interferometers 
in a classical mount, and acquires quasi-monochromatic images 
in the range 560-860 nm (FWHM = 2-4.5 pm). Rapid sequential 
tuning to different wavelengths allows the spectral sampling of a 
number of selected lines, providing spectral information over a 
circular, 80" diameter FOV. Typically, a spectral line is sampled 
in 10-30 spectral points, within a total acquisition time of 3-10 
s. Data can be obtained at the full resolution of 0.083"/ pixel, or 
binned to increase the photon flux and reduce the CCD readout 
time. Coupled with the high-order adaptive optics system of on 
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Fig. 1. Top row: data set 1; bottom row: data set 2. The images refer to data acquired around the middle of the observing sequences. 
Axes are given i n arcsec . Panel s a, e: Broadband continuum at 710 nm. The image in data set 2 has been speckle reconstructed with 
the technique of IWogerl (120071) . Panels b, f: co-temporal HR MDI maps, saturated at + 500 G. Data set 2 attains to an enhanced, 
bi-polar network region. Panels c, g: Intensity in the red wing of the Call 854.2 nm at about 0.1 nm from line core. Bright points 
clearly correspond to small scale magnetic features. Panels d, h: Line core intensity of Call 854.2 nm. Note the large extension of 
fibrils originating in the magnetic elements. 
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Fig. 3. Representative field lines from the MDI-extrapolated potential magnetic field. Left panels: data set 1 . Right panels: data set 
2. The green thick lines trace the open field lines, the red thin line draw field lines closing within the considered area. The white 
dashed circles delineate the IBIS FOV for each data set. The vertical scale in the bottom panels is the height given in Mm. 
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the DST (Rimmele 2004), IBIS images often attain spatial reso- 
lution close to the diffraction limit of the 76-cm telescope. 

2.1. The IBIS data sets 

The data utilized in this work were acquired in two quiet ar- 
eas in close proximity to disk center, on 2004 May 31 (in the 
following data set 1) and 2004 June 02 (data set 2). In both days 
the same acquisition scheme was adopted, sampling sequentially 
the photospheric Fel line at 709.0 nm and the chromospheric 
Call 854.2 nm line. The time necessary to scan the lines was, 
respectively, 4 s and 7 s, while the overall cadence for the full 
sequence was 19 s (the Fell 722.4 nm line was also included 
in the sequence, but not used in this paper). The spatial scale 
was set to 0.166"/pixel (« 120 km at the solar surface). The pho- 
tosphe ric data for set 2 has been utilized in Jansse n~& Cauzzil 
(2006), while both photospheric and chromospheric data, again 
for set 2, have been analyzed in IVecchio et al. d200 7), hereafter 
V07. The Call data for both days has been further described in 
Paper I. 

We analyzed two stretches of 55 minutes of continuous ob- 
servations (175 time steps) obtained in good to excellent seeing 
conditions on each day. Data set 2 had a better seeing than data 
set 1 . We examined the full spectral profiles, as well as studying 
parameters extracted from them. In particular, for each spatial 
position in the FOV we extract by means of spline interpolation 
both the intensity and the position of the minima of the spectral 
line profiles. The latter are interpreted in terms of Doppler-shift, 
with the zero position for the velocity scale defined as the spatio- 
temporal average of the whol e dataset over a quiet po rtion of the 
FOV (see also discussion in iJanssen & Cau zzi 2006, and Paper 
I). 

Fig. Q] gives a synopsis of the data around the middle of 
the observational sequence for both days. The observed FOV 
is shown at different wavelengths together with the simultane- 
ous high resolution (HR) MDI magnetograms. The leftmost col- 
umn shows the broadband continuum at 710 nm, indicating on 
both days a mostly quiet scene, with slightly lower contrast in 
the magnetic regions that delineate the supergranular network. 
At higher magnification these can be resolved as tiny bright 
points embedded within the intergranular lanes (especially in 
the speckle-reconstructed image of data set 2). The second col- 
umn displays the high resolution MDI data obtained simultane- 
ously to the IBIS data sets: the network elements are well dis- 
cernible in the two FOVs, but highlight a very different mag- 
netic environment, with a weaker, unipolar network in data set 
1, and an enhanced bipolar network for data set 2. The third col- 
umn shows IBIS images acquired at about 0. 1 nm from the core 
of the Call 854.2 nm line. Together with the reversed granu- 
lation pattern, these monochromatic images display bright el- 
ements with a one-to- one correspondence to magnetic elements 
dLeenaarts et al.l2006l Paper I). Finally, the fourth column shows 
the images in the core of the Call line, that outline a "segre- 
gated" picture, with part of the FOV occupied by fibrils orig- 
inating from even the smallest magnetic elements, and part of 
it, farthest from the magnetic network, showing an abundance of 
small bright points surrounded by much darker regions (compare 
the temporal evolution of this region from Movie 2 in Paper I). 

2.2. The magnetic environment 

High resolution MDI magnetograms are available for almost the 
whole duration of the IBIS observations on both days. They have 



a spatial scale of 0.6"/pixel, and a cadence of 1 minute. The MDI 
maps outline a very different magnetic configuration of the two 
regions (Fig. Q~]panels b,f). Coronal images confirm this differ- 
ence: while at photospheric level both regions could be classified 
as quiet, EIT Fe XII 19.5 nm images acquired around the time of 
the IBIS observations clarify that the May 3 1 region is localized 
at the edge of an equatorial coronal hole, while the June 2 region 
corresponds to a decaying coronal bright point (Fig. In nei- 
ther case did we observe significant variations of the magnetic 
configuration over the duration of the IBIS observations. 

The differences in magnetic topology are already reflected 
at the chromospheric level, in the morphology shown by the 
Call line core intensities. In set 2 we observe a large number 
of elongated fibrilar structures, with a complex spatio-temporal 
evolution and connecting, in many cases, regions of opposite po- 
larity. On the contrary, the number and length of fibrils in data 
set 1 are much lower, suggesting either short loops that close 
nearby the network elements, or more vertical structures that 
follow magnetic field lines extending into interplanetary space. 
To assess qualitatively the 3-D magnetic field configuration, we 
compute a potential field extrapolation starting from the MDI 
longitudinal field, assumed as the vertical component of the vec- 
tor field at a height of about 200 km. The area utilized for the 
extrapolation is about six times larger than the IBIS FOV. In 
Figure [3] we show the resulting magnetic fields, extrapolated up 
to a height of 2.5 Mm. Representative closed lines within this 
extended region are drawn in red (thin lines), while field lines 
that open to more distant regions or the interplanetary space are 
drawn in green (thick lines). The closed field lines closely re- 
semble the chromospheric morphology: in data set 2 the volume 
is dominated by field lines connecting opposite polarity network 
points, spanning large fraction of the FOV, and reaching heights 
up to about 2 Mm. The quieter data set 1, instead, is character- 
ized by short, and low-lying, closed field lines fanning out from 
the network points and reaching to nearby internetwork weak 
magnetic concentrations. The fibrillar structures so prominent in 
the Ca II core images appear a reliable proxy for the presence of 
a magnetic canopy (Paper I, V07). 

3. Evidence for acoustic shocks in the internetwork 

As for the Call H and K lines, the formation of Call 854.2 nm 
spans a wide range of atmospheric heights (Paper I). From the far 
wings in toward the core, the line gradually samples from the low 
to the high photosphere, while the line core itse lf is formed in 
the low er chromosphere. H owever, as shown bylSkartlien et ail 
( 1994) and more recently bv lPietarila et alj d2006l) . the presence 
of hydrodynamic shocks in the non-magnetic chromosphere de- 
termines a large range of heights for the formation of the line 
core, spanning from as low as 700 km up to 1 300 km. The shocks 
induce rapid variations in the thermodynamics of the affected at- 
mosphere, resulting in asymmetric Ca II 854.2 spectral profiles 
(cf. Figs. 3 and 7 of iPietarila et ai]|2006l) . 

Fig. 2] displays the observed temporal evolution of the IBIS 
Ca II 854.2 nm spectral profiles for four representative pixels, 
two for each data set. The pixels have been chosen as typical of 
the internetwork regions of the FOV. As described in Paper I, 
the sequence of line core images shows a rapidly evolving scene 
with areas that are seething with small bright features immersed 
in an otherwise dark background (compare Movie 2 in Paper 
I). The spectral profiles of Fig. |4]clearly show how this appear- 
ance is the result of the line being repeatedly shifted, strongly 
towards the blue followed by a slower, smaller amplitude drift 
back toward the red. The pattern is repeated at the dominant 
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Fig. 4. Ca II 854.2 nm spectral profiles vs. time (given in s), for four positions within the internetwork. Panels a, b refer to data set 
1; panels c, d to data set 2. The time axis spans the whole duration of the observations. Note the distinctive saw-tooth appearance. 
The thin white line in panel a indicates the evolution with time of the line core Doppler shift. Maximum velocities reach 6-7 km 
s _1 from average position. 



chromospheric periodicity of »120-180 s throughout the whole 
course of the observations. Such strong line shifts are viewed as 
an enhanced brightness when observed in the narrowband, fixed- 
wavelength core images; their spatial coherence determines how 
the bright features are observed with respect to the background 
intensity (see Section[5]l. 

A non-linear steepening of the acoustic waves propagating 
upward from the photosphere at periodicities shorter than the 
acoustic cut-off period of s;180 s is prescribed by the strong de- 
crease in density and the approximate conservation of wave en- 
ergy. This is clearly reflected in the spectra of Fig. [4] where we 
observe both a deviation of the oscillations from a sinusoidal 
form, and a large amplitude of the excursions. At these peri- 
odicities, the typical r.m.s. of chromospheric Doppler shifts is 
around 1 km s" 1 , vs. ~0.07 km s of photospheric motions, cf. 
iReardon et all (2008). However, many of the line shift episodes 
not only display much stronger Doppler amplitudes, of the or- 
der of 5-6 km s _1 , but also exhibit a very abrupt change from 
large redshift to large blueshift, often within a single 19-second 
temporal step in our sampling. Further, the maximum redshift of 
the line often coincides with a brightening episode just blueward 
of the line core, lasting 50-70 s, with enhanced intensities up to 
1 .5-2 times the average value at this wavelength. They are also 
associated with an increased intensity in the extended line wings, 
that progresses in time from the far wings towards the core. For 
the wavelength range of our observations, this time is less than 
about 60 s. Clear examples in Fig. [4] occur at t= 1850 s in panel 
a, or at t=1050, 1200 s in panel d. 

Th ese features are analogous with those of the Ca ll K2V 
grains (lLiull974tlCram et al.lll977tlCram & Darrie1[l98l RU91) 
which have been unequivocally identified as due to hydrody- 
namic shocks in the mid chromosphere (hereafter we will use 
the terms grains and shocks interchangeably). A clear example 



of the resemblance between the temporal evolution of the Ca II K 
and Ca II 854.2 lines in a non-mag netic atmosphere is provide d 
by comparing Fig. [4] with Fig. 2 of Kamio &Kurokawal (l2006): 
all the characteristics previously described for the IBIS data ap- 
pear very obvious in the latter high-resolution observations (note 
the different axes of their Figure with respect to Fig. @}. The 
Doppler shifts are clearer in the Ca II 854.2 line while the bright 
grains and related wing enhancements are less prominent when 
compared to the H and K lines. These differences are a conse- 
quence of the different formation and wavelengths of the two 
lines: the Call 854.2 line has a much stronger Doppler sensi- 
tivity due to its longer wavelength and the Planck function (to 
which the emissivity is related) has a much weaker temperature 
dependence in the near-IR. 



To our knowledge, this is the first time that acoustic shocks 
in the quiet internetwork chromosphere are convincingly (and 
abundantly) observed in the Call 854.2 line, due to both the 
high spectral and spatial resolution of the IBIS data. Indeed, syn- 
thetic profiles calculated with the 1-D hydrodynamical code of 
CS97 well agree with the features descri bed above (Pietarila, pri- 
vate communication; Cauzzi et al. 1 120071) . As further evidence, 
we show in Fig.[5]the development of Ca II 854.2 spectra of two 
example pixels, at several times during the passing of a shock 
wave. The profiles in the right panels are representative of a 
large fraction of internetwork pixels, for which the character- 
istic shifts and asymmetries of the spectral line are in excellent 
ag reement with t h e resu lts of the simulations displayed in Fig. 7 
of iPietarila et al.l ((2006). As an interesting case, we also show 
in the left panels the spectral profiles for a more extreme event, 
where a strong shock greatly deforms the whole line shape up to 
the far blue wing. 
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Fig. 6. Spatial distribution of properties related to shocks. Top row: data set 1; bottom row: data set 2. Results shown in a through 
c, and e through g have been derived from the POD analysis. Axes are in arcsec. Panels a, e: binary "shock maps" at a time around 
the middle of the sequence. Panels b,f. cumulative number of shocks within the whole observational sequence. Panels c, g: maps of 
the shocks' cumulative amplitude over the whole duration of the observations. Panels d, h: spatial distribution of the Fourier power 
for chromospheric velocity, integrated over frequencies at and slightly above the acoustic cut- off (5.5 - 8.0 mHz). 



4. Shock identification 

Typically, shocks in extended FOV have been identified by using 
an intensity threshold on Ca II H^y or K2V filtergrams. This takes 
advantage of the larger temperature sensitivity of the Planck 
function in the blue. However, shocks are foremost a velocity 
phenomena, due to the steepening and rapid reversal of the up- 
ward propagating waves. Hence, the velocity signature is pre- 
sumably the most suitable means to identify their actual occur- 
rence. Given the properties reminded above, Call 854.2 appears 
to be a reliable line with which to identify shocks in the solar 
chromosphere. 

We adopt two different methods that utilize the temporal 
sequence of the full spectral profiles to locate shocks in our 
datasets. For clarity of exposition, we shortly describe here only 
their most relevant characteristics, with further details given in 
the Appendix. The first method we call it the "velocity" method: 
analyzing the temporal profile of the chromospheric line-of-sight 
velocity for each spatial position, we identify times of brusque, 
strong displacements from downward to upward motions. The 
method utilizes a threshold both in amplitude of the displace- 
ment, and in the temporal interval in which this displacement 
occurs, set at a maximum of two time steps of the IBIS sequence 
(i.e. 38 s). The second method uses the "Proper O rthogonal 
Decomposition" technique (POD, Holmes et alj [l996), that sep- 
arates a time-series of spectral profiles into multiple orthogo- 
nal eigenmodes. These are then ordered in terms of decreasing 
energy content, i.e. of their contribution to the whole intensity 
profile. The periodic velocity shifts that produce the sawtooth 
shape in the temporal sequence of internetwork Ca II 854.2 spec- 
tra is clearly identified by the POD in a single antisymmetric 
eigenmode (see panel b of Fig. IA.1I ), properly modulated by a 
temporal coefficient. After finding all the spatial positions for 
which this eigenmode represents the dominant mode, we define 



the shock occurrence at the times when the positive peak of the 
fluctuation shifts from the red to the blue side of the line. We 
apply the same temporal threshold as in the velocity method, but 
no amplitude threshold. 

For both methods we derive a time series of binary "shock 
maps" defining the times of occurrence and spatial distribution 
of the shocks. An example of the POD results is given in Fig.[6]a 
and e. The results from the two methods are similar. As a general 
characteristic, the POD identifies smaller shocks occurring more 
frequently, most notably in the areas filled with fibrils, for which 
the velocity method essentially finds no shocks (see Sects. 15.11 
and [53] below). This can be attributed to the absence of a pre- 
scribed threshold for the amplitude of the shocks in the POD 
case. We will point out other differences and similarities as we 
describe the further results of the analysis. 



5. Shock properties 

5.1. Spatial size and occurrence 

From the shocks maps defined above, we can derive the size of a 
shock by measuring the area covered by contiguous pixels pos- 
sessing values of one. We discarded any shock area covering 
less than 2 pixels in either direction, imposing an effective lower 
limit of 0.3 x 0.3 arcsec 2 to their dimension. For both data sets 
we find average shock extents of 0.4 arcsec 2 (using the POD 
method) to 0.8 arcsec 2 . Under the assumption that the shocks 
are roughly circular, this corresponds to diameters of 0.7 arcsec 
to 1.0 arcsec. These values can be compared to the typical 1-2 
arcsec size of the K2V grains as report ed in the literature (e.g. 
RU91: IOul[l974t iTritschler et"aT]|2007l) . although a number of 
factors set apart our analysis from the more common intensity 
threshold approach. 
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Fig. 5. Temporal evolution of the Ca II 854.2 line during a shock, 
for two internetwork pixels. Intensity is given in arbitrary units. 
The time indicated in each panel is expressed in s from an arbi- 
trary zero point. The thin vertical line indicates the average posi- 
tion of the line core throughout the full spatio-temporal sample. 
Note the broadening of the line during the upward phase. 



The average number of shocks occurring at any time within 
a given area can also be estimated. To this end, we selected for 
both data sets a 30" diameter area roughly covering the center 
of a supergranulation cell centered at (x,y) = (14,0) for data 
set 1 and ((x,y) = (-4,2)) for data set 2 (see Fig. [T}. We find 
an average number of shocks within these areas ranging from 
13 (velocity method, data set 2) to 24 (either method, data set 
1). The number of Kay bright grains in a typical supergranula- 
tion cell have been reported at between 10 and 20 in "the best 
K2V spectroheliograms" (RU91). Our average values are only 
somewhat higher than these previous findings, but we observe 
large variations of the number of shocks during the course of the 
observations, mostly correlated to the seeing conditions: during 
moments of good to excellent seeing, more than 30 shocks can 
be identified within a single cell. 

5.2. Temporal recurrence 

We can further analyze the temporal behavior of the shock oc- 
currence by using two methods. The first one, described in the 
Appendix, shows that the spatial areas where shocks occur in 
large numbers are dominated by shock periodicities around the 
three minutes typical of the chromosphere. We also take an- 
other approach here by calculating the histograms of the "wait- 
ing times" for the shocks, defined simply as the interval between 
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Fig. 7. Distribution of waiting times between shocks, for data set 
1 (panel a) and data set 2 (panel b). 



one shock and the next. No significant difference is found be- 
tween the two methods of shock identification. 

The results are shown in Fig. [7] where a strong peak around 
120-150 seconds is observed for both data sets, with a sharp de- 
cline at shorter time scales and a slower decrease towards longer 
periods. The data do not show any evidence of shock periodici- 
ties below 120 s (i.e. frequencies above 8 mHz), although such 
periods would be observable with our data. The long tail of the 
curves at longer periodicities highlights the fact that Ca II 854.2 
shocks do not occur in a continuous fashion, but often in sepa- 
rate bursts. At the same time, no obvious values of waiting time 
between bursts emerge from the curves, although a change in the 
slope of the distribution appears around 4-5 minutes. Assuming 
this values as discriminant between shocks in the same burst, 
and between different bursts, we find that the average number 
of shocks per burst is only 1.25, and that the average number of 
bursts within the whole sequence is around 3. 

All of this bears on the total number of shocks observed in 
each pixel, displayed in panels b, /of Fig. |6](using the POD re- 
sults). From these maps it is clear that most of the internetwork 
locations only develop shock for a fraction of the observing se- 
quence: only between 25% and 30% of the "shocking pixels" 
display 5 or more shocks during the whole sequence. The maxi- 
mum number of shocks for any given pixel is around 15 for both 
data sets and methods, although this is attained only in a very 
few pixels. 

5.3. Spatial pattern 

A striking result emerging from panels b, fin Fig. [6] is the ex- 
treme inhomogeneity of the spatial distribution of the total num- 
ber of shocks. In particular, large portions of the internetwork 
clustered around the magnetic elements experience very few 
shocks, or none at all. The velocity method provides the same 
spatial patterns, but with even more extended regions displaying 
very few shocks . In the maps of Fig. [6] more than 20% of the 
pixels in data set 1, and 30% in data set 2, do not experience any 
shock throughout the whole course of the observations. These 
values increase to a surprising 50%, respectively 60%, for areas 
with a total number of shocks < 2. These areas are significantly 
larger than the area occupied by the photospheric magnetic ele- 
ments, less than 10% in both cases. These results highlight how 
chromospheric dynamical properties are poorly constrained by 
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Fig. 8. Scatterplot of the cumulative amplitude of shocks vs. 
Fourier power of chromospheric velocities, integrated in a band 
of frequencies just above the acoustic cut-off (5.5- 8.0 mHz). 
Both axes are in arbitrary units. To avoid pixel crowding, the 2- 
D distribution is represented by contours, in logarithmic levels. 



the usually adopted photospheric magnetic diagnostics, a fact 
already remarked in Paper I. 

We further note that these "very few shocks" areas are es- 
sentially coincident with the fibrils observed in the Call line 
core images (compare Fig.|6]with Fig.[TJ, that are again revealed 
as a crucial player in shaping the chromospheric dynamics (cf. 
V07 and Paper I). Conversely, the highest number of shocks is 
recorded in areas well removed from both the photospheric mag- 
netic network and the fibrils. These areas are consolidated in 
small patches scattered within the internetwork. 

5.4. Shocks and Fourier velocity power 

Panels c, g of Fig. [6] provide the cumulative amplitude of the 
shocks over the duration of the observations, again derived from 
the POD analysis (Sect. 1X121 . 

It is interesting to compare the amplitude maps with the dis- 
tribution of the chromospheric velocity Fourier power, integrated 
over the 5.5-8.0 mHz range just above the acoustic cut-off fre- 
quency. The latter is shown in panels d, h of Fig. [6] As discussed 
in V07 for a subset of the data utilized here, in such maps all 
the magnetic elements are surrounded by areas of very low ve- 
locity powe r, corresponding to the "magnetic shadows" first in- 
troduced by Judge et al. (2001) in an analysis of chromospheric 
SUMER data. They are also coincident with the fibrillar struc- 
tures so prominent in Ca II core images and, as clear from Fig. [6] 
with the chromospheric regions where few shocks occur during 
the observational sequence. 

The scatterplots in Fig.|8]display a clear pixel by pixel rela- 
tionship between the velocity power and the cumulative shocks 
amplitude: a simple linear regression finds a correlation of about 
0.65. From these correlations, we surmise that a large fraction 
of what is commonly identified as the "3 -mi nute chromospheric 
oscillations" signal in Fourier analysis (e.g. lOrralll 19661 : [Novesl 
119671: DeubnerjfcFiecJd[l990) is due to the presence of shocks, 
even if they do occur with the intermittent character described 
above. In practice, the linear Fourier analysis identifies the 
shock periodicity and provides a measure of the amplitude of the 
velocity fluctuations at this dominant frequency, even though the 
physical phenomenon that produces this signal is very non lin- 
ear. 
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Fig. 10. Distribution of MDI longitudinal magnetic flux mea- 
sured in points undergoing a large (dotted line) and small (solid) 
number of shocks, respectively. 



6. Shock occurrence vs. magnetic structures 

Numerous observational studies have searched for a connection 
between the occurrence of K21/ grains and the presence of mag- 
netic structures in the internetwork as evidence for a possible 
excita tion mechanism (see e.g. the Introduction of Lite s et all 
1999). The majority of these have concluded that the grains 
are essentially a hydrodynamical phenomenon, for which mag- 
netism plays a minor role at most. However, there remain some 
claims about a tight correspondence between the presence of 
small scal e bipoles in the interne twork and K2V grains, most re- 
cently by Sivarama n et al.l ((2000). Our comprehensive data sets, 
combining high-cadence spectral data over extended areas with 
actual magnetic field measurements, can be used to examine this 
question. 

A coarse inspection of the magnetic maps of Fig.Q]and the 
cumulative shock maps of Fig. [6] immediately shows that shocks 
do not occur at the locations of magnetic structures. A closer in- 
spection shows that the shocks apparently avoid even the small- 
est magnetic elements: in Fig. [9] we show again the cumulative 
shocks maps for both data sets, overlaid with the contours of the 
time-averaged MDI magnetic maps. While in the average MDI 
maps the noise is reduced to a level of about 5 G, in Fig. [9] the 
magnetic contour level is set at 8 G to allow visibility of the 
weak and/or transient structures (the network elements have an 
average magnetic flux of 200-300 G). ^From Fig.|9j one infers 
that even small magnetic structures impede the development of 
the shocks: one notes the small "shadows" in the shock maps 
surrounding several tiny internetwork magnetic structures, for 
example at positions (-30,10); (-10,-7); (11,2) in data set 1, or 
(-14,27); (12,-33) in data set 2. 

A deeper analysis on this issue involves measuring the mag- 
nitude of the magnetic field in those areas that undergo shocks. 
We first selected, in both data sets, the quiet regions where the 
temporal average of the MDI magnetic flux does not exceed 30 
G. We then divided these regions in two parts, one collecting the 
pixels where several shocks (n > 5), over the duration of the 
observations, are observed and the other collecting pixels with a 
smaller number of shocks (1 < n < 4). Excluding the n — pix- 
els effectively avoids most of the magnetic shadows associated 
with the network elements (cf. Fig. [6). The magnetic flux distri- 
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Fig. 9. Cumulative shock maps for data 
set 1 (panel a) and 2 (panel b). Overlaid 
is the contour of the high resolution 
MDI magnetic flux, averaged over the 
course of the IBIS observations. The 
contour level is set at 8 G. Compare 
how even minute magnetic structures 
within the internetwork, probably oc- 
curring only for a fraction of the ob- 
servations, correspond to a decreased 
number of Ca II shocks. Spatial scale in 
arcsec. 



butions obtained from the temporally resolved MDI maps in the 
two classes of pixels are shown in Fig. [10] The vertical line at 15 
G represents the noise level for a single pixel in MDI maps. In 
both data sets, the distributions for the two classes of pixels start 
diverging shortly after the noise value, with the pixels undergo- 
ing many shocks characterized by sensibly lower values of the 
magnetic flux. 

The curves diverge at the level of 1CT 2 in relative occur- 
rence, i.e. only a small fraction of the pixels displays this effect. 
Nevertheless, the difference is clear. Combined with the maps 
in Fig. [9] this shows that the internetwork magnetic fields actu- 
ally seem to have a negative effect on the development of shocks. 
Magnetograms with higher sensitivity and spatial resolution than 
those of MDI might reveal if there is a critical area or flux value 
below which this effect no longer holds. 

7. Shock occurrence vs. photospheric velocities 

We now examine the relationship between the occurrence of 
shocks and the underlying photospheric dynamics. We do so 
both with classical Fourier analysis, and also using wavelet anal- 
ysis on photospheric velocities, which is a more efficient way to 
take into account the highly intermittent character of the shocks. 

7.1. Fourier phase difference spectra 

The spatially-averaged Fouri er velocity power sp ectra for these 
data sets are displayed in iReardon et al.l d2008l) and Paper I, 
and are in general very similar to previously published spec- 
tra (e.g. Deubner & Fleck 1990): the photospheric curve shows 
a strong maximum around the classical 5 minute periodicity 
(v ~ 3.5 mHz), while the chromospheric curve for the internet- 
work shows a flatter maximum from ~ 4.5 to 7 mHz. Here we 
concentrate instead on the phase relationship between the pho- 
tospheric and chromospheric velocities as a means to study the 
propagation of acoustic waves. 

Panels a, b in Fig. [TT] display the phase difference spectra 
between the Fourier transforms of the chromospheric and photo- 
spheric velocities, obtained over the pixels that develop at least 
5 shocks throughout the observational sequence. The Figure is 
built by plotting the binned phase difference A(f> weighted by the 
cross-power amplit ude and normalized per temporal frequency 
bin; for details see Krii ger et al.l d2001l) and references therein. 
The noisier scatterplot for data set 1 is due to the worse seeing 
conditions with respect to data set 2. A positive phase difference 
describes signals propagating upward in the solar atmosphere: 



the scatterplots thus clearly evidence the propagation of acous- 
tic waves from photosphere to chromosphere, up to about 10-12 
mHz. The bottom panels c, d provide the coherence spectra of 
the Fourier velocity signal (solid line). The coherence is very 
high around the velocity power peak at 5-6 mHz (see Paper I), 
and falls below the confidence limit of 0.5 around 8-9 mHz. 

If we consider an equal-size sample of pixels characterized 
by a smaller number of shocks, we obtain essentially the same 
phase differences (not shown in Figure) and slightly lower coher- 
ence values (dashed line): the same photospheric piston mech- 
anism thus seems at the base of all quiet chromospheric os- 
cillations, regardless of whether other atmospheric properties 
are conducive to the development of many mid-chromospheric 
shocks or not. The coherence is instead much smaller if we use 
an equal sample of pixels belonging to the fibrillar areas (dotted 
line); these are areas where most pixels do not develop any chro- 
mospheric shocks, especially for data set 2. In these atmospheric 
regions, a direct (vertical) relationship between photospheric and 
chromospheric dynamics no longer holds. 

The coherence spectra display a curious bump around 7 
mHz, most evide nt in data set 2. S uch a feature has already 
been reported by Deubner & Fleck (1990) in a spectrographic 
analysis of a quiet region. Like these authors, we find no obvi- 
ous explanation for the feature. We suspect that it might be re- 
lated to the "aureoles" of enhanced high frequency power some- 
times measured in both photospheric and chromospheric signa- 
tures around active regions or strong network elements (see e.g. 
iKrijger et al.ll200U and references therein). The bump is indeed 
more evident in data set 2, which has stronger magnetic fields 
and a more "active region-like" magnetic configuration than data 
set 1. 

7.2. Photospheric velocity power and shock occurrence 

In their analysis of K2V grain formation, CS97 remarked that 
photospheric acoustic waves at or near the cut-off frequency play 
the most important role for the development of shocks. In par- 
ticular, they stated that whenever photospheric velocities have 
significant power around 5 mHz, there will usually be signifi- 
cant bright grains. To check these statements, we performed a 
wavelet analysis of the photospheric velocities. We concentrate 
on the spatio-temporal features showing enhanced power at pe- 
riodicities between 120 and 200 s (5-8 mHz), taken as values 
larger than the average over the whole sample (this accounts for 
about 35% of the total number of pixels). Hereinafter we refer to 
them as "areas", even if they contain a temporal dimension. 
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Fig. 11. Panels a, b: Fourier spectra of phase differences between 
the photospheric and chromospheric velocities, for data set 1 and 
2, respectively. The diagram has been calculated over the points 
in the FOV that develop at least 5 shocks throughout the obser- 
vational sequence. Panels c,d: corresponding coherence spectra 
(solid line). Dashed lines refer to the coherence measured for 
an equally populated sample of pixels undergoing less than 5 
shocks; dotted lines to the coherence for the case of fibrils. 



Fig.fl2lshows example that the mechanism outlined by CS97 
is most probably at work. In all panels the x-axis represents the 
spatial dimension along a horizontal cut within a quiet region in 
the FOV, while the y-axis displays the temporal dimension over 
the duration of the observations. Top panels refer to instances 
in data set 1; bottom panels to data set 2. The small diamonds 
represent the occurrence of Ca II 854.2 shocks, while the con- 
tours outline the areas of enhanced photospheric velocity power 
(left and center panels), and the presence of magnetic features 
(right panels). In the left panels, most of the patches of strong 
photospheric velocity power encompass the areas where shocks 
are observed, indicating a relationship between these two phe- 
nomena. The coincidence in the case of data set 1 is particu- 
larly striking, with photospheric contours enclosing the chromo- 
spheric shocks for the whole sequence over extended areas. We 
note that in order to obtain this coincidence a delay of about 
120 seconds (6 time steps) has been applied between the photo- 
spheric signal and the occurrence of shocks. This is consistent 
with the idea that the photospheric perturbations propagate up- 
ward at the speed of sound over a ~ 700-1000 km height differ- 
ence. 

However, determining whether all areas of enhanced photo- 
spheric power indeed give rise to shocks, proved a difficult task. 
As the presence of magnetic fields breaks the direct relationship 
between the photospheric and chromospheric dynamics, one has 
to carefully select the areas to scrutinize this phenomenon. The 
center panels of Fig. [12] provide a clear example of this. The 
plots have been obtained for spatial positions still well removed 
from the magnetic network, and in them we observe that in sev- 
eral instances, especially for the data set 2, areas of enhanced 
photospheric power do not produce any shocks, contrary to the 
expectations of the hydrodynamical modeling. The small mag- 
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Fig. 13. Distribution of photospheric velocity power between 
120 and 200 s periodicity for positions undergoing shocks (thick 
curves) and not (thin curves). Panel a refers to data set 1; panel 
b to data set 2. The overall higher power values in data set 2 are 
due to better seeing conditions. 



netic structures responsible for this lack of shocks are outlined 
in the right panels, at a level of ~ 30 G: it is obvious how the 
shocks generally avoid both the magnetic elements proper, and 
some surface around them, as for the larger network shadows. 

We hence attempted to identify the quietest portions of the 
FOV for both data sets, by combining the series of MDI maps 
with the Call 854.2 core intensity maps (see also Fig. I141>. In 
both cases we selected a 8"xl3" area (50x80 spatial pixels, con- 
sidered independent) positioned almost at the center of the FOV, 
that displayed the least magnetic flux in MDI maps and no ob- 
vious evidence of fibrils in Call core images. For these regions, 
we find that more than 80% of the areas with enhanced pho- 
tospheric velocity power contain shocks. Further, the number 
of shocks measured in each area is directly proportional to the 
spatio-temporal extent of that area: the development of acoustic 
shocks appears to be continuos as long as photospheric condi- 
tions are favorable. As a rough estimate of the influence of mag- 
netic structures, we compare the total number of shocks mea- 
sured in these quietest portions of the FOV, taken as prototypes 
for the undisturbed acoustic shocks mechanism, with those mea- 
sured in the whole supergranular cells defined in Sect. 15.11 We 
find that the quietest regions produce, per unit area, from almost 
2 (data set 1) up to 3 (data set 2) times more shocks than the 
supergranular cells as a whole. 

We have shown that high photospheric velocity power at the 
cut-off periodicity is certainly responsible for the development 
of shocks. But do all shocks occur in areas of high power? To 
answer this question, we show in Fig.[13]fhe distribution of the 
photospheric velocity power for both data sets. The curves have 
been calculated for all the pixels undergoing shocks (thick line) 
and pixels where shocks are not observed (thin line). Again a 
delay of 120 seconds has been applied between photospheric 
and chromospheric signals. The photospheric power is obviously 
stronger in the former case than in the latter, as expected from 
the results described above, with an average value 45% larger 
for data set 1, and 35% for data set 2. By measuring the area en- 
closed below the curves, we find that in both data sets over 55% 
of shocks occur in areas with photospheric power larger than the 
average value (the latter accounting for 35% of the total pixels, 
as stated before). This is to be considered a lower limit, as shocks 
that are part of a common pattern sometimes lie just outside of 
the areas defined by the photospheric level, as clearly seen in Fig. 
[T2l Thus, a large fraction of shocks is due to enhanced power in 
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Fig. 12. Examples of spatio-temporal occurrence of chromospheric shocks (small diamonds) in relation to the photospheric velocity 
power at periodicities between 120 and 200 s (left and center panels) and magnetic features from MDI (right panels). Velocity power 
contour levels are set at the average value over the whole FOV plus half standard deviation. Top row refers to data set 1 ; bottom row 
to data set 2. The temporal axis is shorter than the full IBIS observational sequence as wavelet maps are not reliable in the initial and 
final temporal samples, and because of the delay applied between photospheric and chromospheric signal. The x - axes coordinates 
are given in the same units as in the FOV of Figs. Q] [6] Left panels correspond to y— 5, for data set 1 and data set 2, and define 
quiet areas with strong correspondence between photospheric dynamics and shocks' occurrence. Center panels correspond to y= 12, 
3 for data set 1 and data set 2, and show areas with comparable photospheric dynamics but notably reduced development of shocks. 
This is due to the presence of small-scale magnetic structures, as outlined by the contours in the right panels. Contours are set at 30 
G for the magnetic flux. For data set 1 the MDI sequence ends before the respective IBIS temporal sequence (cf. upper right panel). 



the photosphere at periodicities near the cut-off. Possible causes 
for shocks not accounted by this mechanism could be the addi- 
tional contribution from high-frequency waves when waves near 
the acoustic cutoff frequency are weak (CS97), as well as some 
instances of non-vertical propagation. 



8. Summary and Discussion 

The temporal sequences of internetwork Ca II 854.2 nm spectra 
show compelling evidence for abundant acoustic shocks in the 
mid-chromosphere (Sect. |3J. The evolution of spectral charac- 
teristics of the Call 854.2 line during these events is in com- 
plete analogy with the more famous Call K.2V grains, clearly 
explained by CS97 as due to weak acoustic shocks in the mid- 
chromosphere. We have set out to provide a comprehensive view 
of the occurrence of these shocks, by combining a large ar- 
ray of complementary diagnostics, including spectral informa- 
tion in both photospheric and chromospheric lines; high spatial 
resolution over an extended field-of-view; simultaneous high- 
resolution magnetic maps from MDI; and observations of two 
different quiet Sun targets with different magnetic topologies, 
for a period of about 1 hr each, at high cadence. We surmise 



that such a dataset is rather unique, and indeed provided much 
insight into the issue. 

Shock identification. We have identified the spatio-temporal lo- 
cations of chromospheric shocks by using a velocity criterion on 
the Call 854.2 nm line, namely the abrupt displacement of the 
profile from red to blue. This approach provides several intrin- 
sic advantages ov er the more commonly used intensity threshold 
methods (see e.g. lTritschler et ai1 l2007, and references therein): 
Using a velocity determination avoids confusion between bright- 
enings from shocks and those due to transient small-scale mag- 
netic structures within the internetwork. The latter correspond 
to what has bee n termed "magneti c grains" (RU91) or "per- 
sistent flashers" dBrandt et alJll992t Ide Wiin et al.ll2008l) . and, 
while as bright as the internetwork grains, they are character- 
ized by longer evolutionary timescales. Analysis of filtergrams 
in Ca II H and K may confuse the intensity signature from these 
two different processes. Our shock identification technique also 
identifies the start time of each shock, assuring that they are 
counted only once during the temporal sequence, something that 
can be difficult with threshold methods. Finally, by using the 
full spectral profile, we reduce the influence of spatially and 
spectrally scattered light and photometric errors. The velocity 
method for identifying shock relies on a threshold for the am- 
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plitude of the red-to-blue Doppler shift, but while some of the 
shock parameters did show a dependence on the threshold value 
(in particular the size, and total number), others, such as the tem- 
poral characteristics and spatial distribution, did not. 

Shocks' properties. Many of the properties that we derive for 
the Call 854.2 shocks are analogous to those reported in the 
literature for K.2V grains, confirming that the two phenomena 
represent the same physical process. In particular, we find that 
the Call 854.2 shocks often appear in bursts, with the typical 
interval between shocks ranging from ~2 to 4 minutes (peak- 
ing around 150 s), and intervals between bursts of 5 minutes 
and longer. This intermittent behavior is a reflection of their es- 
sentially stochastic nature (CS97). No obvious periodicity below 
120 s is observed. 

The shocks are not ubiquitous, neither in space nor time. 
Only 25 to 30% of the quiet internetwork displays five or 
more shocks during the ~ 1 hour sequences. Assuming that 
shocks giv e rise to K.2V em ission for about 60-80 s for each 
event (e.g. iBeck et al.ll20 08). this rate of occurrence translates 
into a measurable K~2y emission for about 6-10% of the total 
spatio-temporal sample, in agreement with earli er works (e.g. 
Ivon Uexkuell & Kneerlll995l:(steffens etaTII 19961) . Recently, us- 
ing high resolut ion Call H spectra acquired with POLIS, 
IBeck et ail d2008h claimed that in quiet Sun regions the core of 
the Ca II H resonance line spends about half of the time in emis- 
sion. This is probably misleading, as they implicitly assumed 
that shocks occur everywhere and continuously with the typical 
3 minutes periodicity, an assumption obviously not fulfilled on 
the Sun. 

The average shock size of 0.7"-1.0" arcsec measured from 
the Call 854.2 data is sensibly smaller than the l"-2" arcsec 
typically reported for K.2V grains. We believe our measurement 
comes closer to the true size of the shocks due to the higher spa- 
tial resolution reached with IBIS with respect to older works, 
the limited influence of instrumental scattered light on our mea- 
surements, and the reduced scattering in the solar atmosphere 
for the Ca II 854.2 line with respect to the H and K. A small ex- 
tension of the shocks explains why earlier spectral observations 
of the Ca II 854.2 nm in the quiet internetwork, obtained with a 
0.6" pixel size, did not display the clear sawtooth behavior such 
as that seen in Fig.|4](cf. Fig. 3 in lDeubner & Flecklll99dl) . The 
smaller effective size of the Ca II 854.2 shocks, together with the 
overall better resolution, accounts for the larger occurrence rate 
that we measure in a typical supergranular cell (>20 in our case, 
vs. 10-20 in the "best spectroheliograms", cf. RU91). 

Spatial distribution and magnetic shadows. An important re- 
sult of our analysis is the extreme inhomogeneity of the spatial 
distribution of the shocks (Sect. |53] l. It has long been known that 
shocks can be relati vely rare and wi dely dispersed in the super- 
granular cells (e.g. Lite s et ail 1 19941) . Indeed we find that up to 
50-60% of the FOV display zero or very few shocks, even for the 
case of data set 1 which, being situated at the edge of a coronal 
hole, arguably has one of the weakest possible magnetic config- 
uration. However, with our data we clearly show that the shocks 
concentrate in the internetwork but away from the regions where 
chromospheric fibrils are present. The local magnetic topology 
outlined by the fibrils is very different for the two data sets, with 
many short loops closing back from the network elements to the 
near internetwork regions in set 1 and longer and more stable fib- 
rils nearly spanning the supergranular cell in set 2 (Fig. [3}. This 
nature of the canopy and the spatial distribution of the shocks is 
not easily predicted using only photospheric indicators such as 



the granulation in the continuum, or the location of the network 
magnetic elements. 

A second interesting result is that the total amplitude of 
the shocks over the course of the observations is directly re- 
lated to the magnitude of the Fourier power of chromospheric 
velocities with periods around 3 minutes (Fig. [8}. Thus, the 
shocks represent an important, if not dominant, c ontribution to 
the "3- minu te chromospheric oscillations" (e.g. IOrrall [l966: 
lNoveslfl967h . even if they do occur with an intermittent char- 
acter. Because the shocks are related to large excursions of the 
spectral line, corresponding to several km s when parametrized 
in terms of line core Doppler shifts, they can contribute a large 
fraction of the velocity signal around the acoustic cutoff fre- 
quency in classical Fourier analysis, even if they do not occur 
continuously. The two quantities are essentially equivalent, at 
least in the quiet Sun. 

These two results fit in nicely with several pieces of evi- 
dence reported in the literature about the so-called "magnetic 
shadows", and allow us to for m a coheren t expla nation of their 
existence. First introduced by Judge et al. (2001), the shadows 
are regions of reduced Fourier oscillatory power at the 3-minute 
periodicity, often observed in low- and mid-chromospheric sig- 
natures around (but not coincident with) magnetic network el- 
ements. They have been reported in intensity obse rvations of 
the c ontinuum at 119 nm obtained with SUMER ( Judg e~et~ai1 
2001 ), and of the c ontinuum at 160 nm obtained with TRACE 
dKrijger et al.ll200"Tb . as well as in Fourier anal ysis of chromo- 
spheric velocities by V07 and lLites et alj d 19941) . 

All these works invoke the presence of horizontal magnetic 
field lines (the "canopy") to somehow disturb the propagation 
of normal acoustic waves towards the outer solar atmosphere, 
for example by means of wave scattering or mode conver- 
sion ([Judge et alj2001l:rMcIntosh & Judgel2001l:[McIntosh et al.l 
2003). With our data, we prove that the shadows correspond to 
a strong suppression of the number of acoustic shocks. Further, 
we show that the shadows are tied to the presence of chromo- 
spheric fibrils, i.e. to structures with a very distinct physical 
structure with respect to the surrounding atmosphere. Indeed, 
we believe they are best compared to the Ha mottles and fib- 
rils observed around quiet Sun magnetic network, that can be 
explained either by elevated structures s uch as spicules , or "em- 
bedded" ones (see e.g. the discussion in Al et al. 2004). Hence, 
the upward propagating acoustic waves may be disturbed not so 
much by the horizontal field itself, but by the different stratifica- 
tion of the atmosphere. This should be taken into account in the 
discussions in terms of high- or low- plasma /3 (the ratio of ki- 
netic to magnetic pressure) as the main discriminant between an 
'acoustically-dominated' and a 'magnetically-dominated' atmo- 
sphere. An estimate of /3 that uses the classical, semi-empirical 
temperature and density values (e.g. from VAL) might not be 
appropriate in these regions. In the case of elevated structures, 
it is also possible that the acoustic shocks develop normally in 
the low chromosphere whenever photospheric conditions are fa- 
vorable, and the higher-lying fibrils simply act as a mask, hiding 
the lower laye rs from view. Some support for this scenario can 
be gathered in lde Wiin et alj d2007l) as well as from Movie 2 in 
Paper I, where the lateral motions of the fibrils at times "un- 
cover" the presence of bright grains, especially towards the end 
of their horizontal extension (see also Fig.[T4b. 



Photospheric drivers. The analysis of Sec. 17.11 clearly shows 
that the Ca II 854.2 shocks partake in the general chromospheric 
dynamics, responding to acoustic waves propagating from lower 
layers. The phase relationships between photospheric and chro- 
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Fig. 14. Maximum (top row) and minimum (bottom row) intensity of the Call line core throughout the 55 minute sequences. 
Left: data set 1; right: data set 2. The scale has been set in radiation temperature (K, see text for details). This representation of 
the data immediately provides a distinction among the three major chromospheric components: the bright network; the acoustic 
chromosphere where acoustic shocks show up as tiny bright points; and the magnetic chromosphere threaded by slanted fibrils that 
can appear as either dark or bright. 



mospheric velocities in regions undergoing shocks are fully con- 
sistent with earlier results obtained for general internetwork ar- 
eas, proving the propaga tion of acoustic waves up to 10-12 mHz 
dDeubner & Fleckl fl990). The occurrence of shocks in this gen- 
eral oscillation pattern must then be triggered by complementary 
circumstances. By means of a wavelet analysis we find that in- 
stances of high photospheric velocity power at periodicities be- 
tween 120 and 200 s (5-8 mHz) are well correlated with the 
occurrence of chromospheric shocks. Tellingly, the correlation 
is maximized when allowing for a delay of ~120 s between 
the photospheric and chromospheric signal, indicating a photo- 
spheric perturbations propagating upward at the speed of sound 
over a ~ 700-1000 km height difference. 

The small dimension of the shocks is not completely corre- 
lated with the piston's size - the analysis of Sect. |7.2| shows how 
areas of coherent photospheric oscillations at the relevant peri- 
odicity have lateral sizes of up to several arcsec. However, it is 
important to note that these areas completely encircle the ma- 
jority of the pixels that develop shocks (left panels in Fig. [12} , 
All of these results are consistent with the analysis of CS97 that 
state how the appearance of the bright grains, including their 
size, is fully determined by powerful waves near or just above 
the acoustic cutoff frequency interfering with higher frequency 



waves. The interference pattern that determines the actual size 
and recurrence of shocks will be contained within, but not nec- 
essarily coincident with, the regions of enhanced photospheric 
velocity at the dominant frequencies. 

The strength of photospheric motions is not the only fac- 
tor shaping the "quiet" chromosphere and its dynamics. We find 
that the presence of magnetic structures has a profound influ- 
ence on what occurs in the lower chromosphere and, in particu- 
lar, whether the acoustic mechanism outlined above can operate 
undisturbed. Many areas with suitable photospheric dynamics 
do not develop shocks (the "shadows") and show a decreased 
coherence between the vertical photospheric and chromospheric 
velocity signals at all frequencies (panels c and d of Fig. ITTb . 
The dynamics of the fibrils present in the shadows' areas prob- 
ably is most-closely related to the oscillatory behavior at their 
footpoints than to the directly underlaying layers. Interestingly, 
we note that pixels developing a small number of shocks have 
a lower coherence than the ones developing many shocks; we 
assume that this is related to a mechanism similar to the shad- 
ows, just operating on smaller scales. Several hints support this 
hypothesis, such as: the diffuse appearance of the Ca II line core 
images (Fig. Q] especially panel d) that might betray short, un- 
resolved fibrils; the correlation between the presence of small, 
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transient magnetic structures observed in the MDI HR maps and 
the reduction in total shocks numbers (Pig. [9j»; the lack of chro- 
mospheric shocks in areas corresponding to enhanced photo- 
spheric velocity power patches, if there exist nearby small scale 
magnetic structures (Sect. lTTZl l. 

Vertical propagation. Throughout this paper, we searched for 
correlations between photospheric and chromospheric signal un- 
der the assumption of vertical propagation. This is of course a 
limitation of our analysis and has been invoked often as one 
of the major problems i n the CS97 study of K2V grains (e.g. 
lUlmschneider et aTl 120051) . However, there are several indica- 
tions that this assumption is not misguided. First, the level of 
coherence between the photospheric and chromospheric veloci- 
ties at the dominant frequency of ~5 mHz is very high, especially 
for the quieter data set 1 (Fig.fTTTi. This would not be the case if 
the propagation were not predominantly vertical. Indeed, we find 
a much lower coherence value for the case of the chromospheric 
fibrils, which have an obviously slanted geometry. Second, the 
wavelet analysis of Sect. 17.21 shows how patches of high photo- 
spheric power at the relevant periodicity well coincide with the 
appearance of shocks in the chromospheric areas directly above. 
As the photospheric patches of high power extend for 1500- 
3000 km in the horizontal dimension, i.e. at least twice as much 
as the vertical separation between the regions of formation of 
the Fel and Call lines, a small deviation from the vertical (say 
below 30°) would also not cause any drastic difference in our 
results. 

The wavelet analysis shows that in some instances chromo- 
spheric shocks are measured even in absence of high photo- 
spheric power; these could be obvious candidates for a search 
on the effects of inclined propagation, and deserve further in- 
vestigation. They could also be related to the presence of high 
frequency waves or photospheric power just below our cutoff for 
enhanced power. Still, it appears that small magnetic structures 
represent a more important factor in the shaping of the 3-D chro- 
mospheric dynamics, reducing the overall number of shocks by 
large amounts (Sect. 17721 1. 

Three quiet chromospheres. From our results, we propose that 
the "quiet" lower-chromosphere should be considered to com- 
prise at least thre e components, as somewhat anticipated in 
Judg e et afj d2001l) and remarked in Paper I. These three com- 
ponents are: 1- the purely acoustic chromosphere, where the 
strong photospheric motions at frequencies near the acoustic cut- 
off can propagate upwards undisturbed and develop into shocks; 
2- the network and internetwork magnetic elements, that open 
up partially from the photos phere and are heated by a m echa- 
nism still unexplained (see lHasan & van Ballegoo iien 2008] and 
references therein); and 3- the inclined fibrils that magnetically 
connect horizontally separated areas. The fibrils can span as long 
as a supergranular cell (as in data set 2) or as short as few arcsec- 
onds (as in the smallest shadows in Fig.[6](f); in any case they de- 
fine chromospheric volumes with very different dynamical prop- 
erties than the other two components. In particular, they show a 
reduced chromospheric vertical velocity power at all frequen- 
cies (see also Paper I); distinctly longer evolutionary timescales 
( V07); a lack of vertical coherence between the photospheric and 
chromospheric signals (Sect. 17711 1. 

These three components can be immediately identified by 
using maps of the Call 854.2 line core intensity, a reliable indi- 
cator of temperature (Cauzzi et al., in preparation). In particular, 
we show in Fig.[14]the maximum (top row) and minimum (bot- 
tom row) value of the line core intensity, measured for each spa- 
tial pixel within the 55 min observing sequences. In the Figure 



the intensity has been converted to radiation temperature by scal- 
ing the average intensity computed over t he full spatio -temporal 
sample to the calibrated atlas intensity of lNeckell d 19991) and ap- 
plying the inverse Planck function (the average intensity corre- 
sponds to 3800 K). Obviously this definition does not represent 
the real kinetic temperatures as the line is not formed in LTE, but 
it does provide an indication for the magnitude of the variations. 
In Fig.[14]the acoustic shocks stand out clearly in the maximum 
intensity maps as bright specks with a quasi-granular appearance 
(but with smaller dimensions) and temperatures up to 4200 K. At 
times, especially in data set 2, the shocks can even be discerned 
in amongst the long fibrils, supporting the idea of the latter op- 
erating as a switch in between the lower chromosphere and the 
observer. The canopy-free internetwork regions where there are 
significant numbers of shocks are also the regions of the low- 
est intensity, and presumably lowest temperature, in the entire 
sequence, with radiation temperature values as low as 3400 K. 
The magnetic elements appear also as distinct point-like sources, 
but larger and brighter than the acoustic shocks, with tempera- 
ture excursions ranging from a minimum of 4000 to maxima of 
4700 K. Finally, the slanted fibrils have a very different morphol- 
ogy and might be better defined as the "diffuse" sources in both 
the maximum and minimum intensity maps. Two things about 
them are particularly remarkable: they are the component whose 
intensity (temperature) changes the least between maximum and 
minimum values, with largest values close to the network points 
and cooler ends towards the internetwork; and they cover a large 
fraction of the FOV, even for data set 1 . The magnetic influence 
appears very pervasive over the quiet Sun. 

Finally, it would appear that Fig. [14] also offers a defini- 
tive statement on the issue of the spatial correspondence be- 
tween small scal e magnetic structures and Y^v bright points (e.g. 
ISivaraman et al. 2000). Indeed, many such small magnetic con- 
centrations are visible in the internetwork portion of the FOV in 
both of our data sets, and they do correspond to enhanced bright 
"points", some of them only 1-2" wide. However, they also are 
clearly distinguished from the acoustic grains by the diffuse co- 
coon around them, made up most probably of unresolved fibrils. 
They a ppear just like min iature versions of normal network ele- 
ments dKriiger et al.l200lb . The analysis of Sect.|6]makes it clear 
that these bright structures do not participate in the typical dy- 
namics of acoustic shocks, making them K2V bright "points", but 
not K.2V bright grains! Claims to the contrary, such as those of 
ISivaraman et al.l d2000l) . stem from the limited coverage of the 
available datasets or the limited number of co-spatial magnetic 
and K2V measurements. 

9. Conclusions 

In our analysis, we reach two main conclusions, both of which 
are represented in Fig. [14] 

The first is that temporally resolved spectra of the Ca II 854.2 
nm provide a comprehensive picture of the occurrence of acous- 
tic shocks in the quiet solar atmosphere. We determine that these 
shocks are a fundamental part of the general chromospheric dy- 
namics, as they essentially make up the bulk of the three-minute 
oscillations in the middle chromosphere. We find that the shocks 
occur in direct response (after a 120 second delay) to power- 
ful photospheric oscillations near the acoustic cutoff frequency, 
with a primarily vertical propagation. No major role is found for 
waves with periods below 120 seconds or with a slanted prop- 
agation, c onsistent with the 1-D p urely hydrodynamical simu- 
lations of ICarlsson & Stein] dl997l) . Portions of the quiet chro- 
mosphere are then the site of brief episodes when the acoustic 
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shocks dramatically enhance the temperature (Fig. [14] a and b) 
mixed with extended periods of lower temperatures (Fig. [141c 
and d), as opposed to the steady-state mean temperature rise typ- 
ically assumed in semi-empirical mod els (see the discussion in 
ICarlsson & Steidll995tlCarlssonll2007l) . 

The second, more important conclusion is that these acoustic 
processes are very often significantly disturbed by the presence 
of magnetic field, even in the quiet Sun. Both at the network and 
(surprisingly) at internetwork scales, the signature of acoustic 
shocks is suppressed both in and around the footpoints of mag- 
netic structures. Hence the distribution of shocks is considerably 
reduced compared to what would be inferred from the amplitude 
of the photospheric oscillations alone. The main player in this 
respect appears to be the magnetic topology, and in particular 
whether the field has an open or closed configuration. The ar- 
eas of modified acoustic behavior coincide with the presence of 
chromospheric fibrils, betraying an atmosphere that is stratified 
differently with respect to the classical 1-D view. 

These results provide an important hint concerning the 
longstanding controversy about why the dynamical model of 
CS97 fail s to predict the ubiquitous emission observed in UV 
lines (e.g. Kalkofen etaill 19991: lKalkoferil200lt ICarlssoril2007t 
lAvrett & LoeserH2008l) . Our analysis is based on observations in 
the Call 854.2 nm line, a spectral signature formed relatively 
low in the atmosphere. Our data provide uncontroversial evi- 
dence that even at these heights the quiet atmosphere is highly 
structured by pervasive magnetic fields that connect different lo- 
cations. The presence of fibrils is a clear testimony of these inter- 
connectiono It is then natural to assume that other spectral di- 
agnostics forming in higher, less dense parts of the atmosphere, 
experience such a ma gnetic structu ring to an even larger extent, 
much as advocated bv lAvrei d2002l) . 

Far UV lines and continua are one such diagnostic, and they 
figure prominently in the derivation of the temperature structure 
in se mi-empirical models, o ften as spatial and temporal averages 
(e.g. lAvrett & L oeser 2008). It follows that a substantial fraction 
of the UV flux could originate in magnetic structures, even in the 
quiet internetwork chromosphere. This could happen either be- 
cause of the volume filling of the magnetic field at these heights, 
or because of the dominant weight of magnetic -related bright- 
ness in the averages: Fig.[l4]illustrates how magnetic concentra- 
tions introduce a brightness component with a different spatial 
distribution and temporal evolution than for the purely acoustic 
case. This happens on spatial scales at the limit of the resolution 
of many modern instruments, including SUMER. Any averaged 
or marginally resolved observation of the chromospheric flux 
might then contain a sizable magnetic component, in ways not 
immediately predictable by using only photospheric indicators. 
A dominance of the magnetic topology in the quiet upper chro- 
mosphere would also explain why UV observations obtained 
with SUMER have produced so many disparate results. These 
range from the cle ar signature of upw ard propagating acoustic 
waves reported by Wik st0l et al.l d2000t) all the way from the 
photosphere up to the base of the corona, as could happen in the 
case of a coronal hole, to instances of complete disconnection 
between the observed low er and upper chromospheric dynamics 
(Judge etaUbOOll 120031) . as would be the case for regions with 
a locally closed magnetic topology. In this scenario, obviously 
there is no need to reconcile the predictions of CS97, a purely 



1 We note that this magnetic structuring has been observed for 
decades in Ha, a diagnostics that has figured little in the discussion 
on chromospheric heating. 



non-magnetic model, with the UV observations, as they pertain 
to completely different physical regimes. 

We thus conclude that the radiative losses in the "non- 
magnetic" (quotes are necessary at this point) solar chro- 
mosphere are strongly influenced, and possibly dominated, 
by processes related to magnetic fields. We remark that this 
same conclusion has been reached by Judge and collaborators 
(lJudge & Carpenter|[T998l:lJudge et atll2003l 120041) , in a precise 
analysis of UV emission lines for both the case of the Sun and 
other stars with convective envelopes. The presence of fibrils 
in large fractions of the quiet atmosphere provides a plausible 
agent for such dominance. On the one hand, they disrupt the 
normal vertical propagation of acoustic waves from lower lay- 
ers, strongly reducing the effects of a purely acoustic mech- 
anism in shaping the chromosphere. On the other, they pro- 
vide the means to transport energy of magnetic origin, or me- 
diated by the magnetic elements, from the network footpoints 
towards the center of the supergranular cells. We note that this 
idea is consistent with several recent works that have called 
into question whether the photospheric acoustic flux at high fre- 
quencies i s sufficient to compensate the c hromospheric radia- 
tive losses dFossum & Carlssonll2005[ [20061) or, even more perti- 
nent, whether different energy sources such as low- frequencies 
magneto-acoustic waves, or atmospheric gravity waves, might 
play a more important role ( Je fferies et al.l 12006c IStraus et all. 
120081) . 
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Appendix A: Identification of shocks 

A. 1 . The velocity threshold approach 

A way to identify acoustic shocks "events" is through the sawtooth peaks in 
the velocity field, defined via the Doppler shift of the line intensity minimum. 
The small white crosses in Fig. [4] panel a show this parameter for the particu- 
lar pixel. As apparent from both Fig.|4]and|5] during the upward phases of the 
shock the spectral line broadens considerably in the blue wing, probably lead- 
ing to an underestimate of the actual value of the vertical velocity by using the 
simple Doppler shift of the minimum. However, such an effect does not seem 
particularly severe (compare for example Fig. 7 in Pietarila et al. 2006), and the 
sign of the recovered velocity is always correct. 

Following the temporal behavior of the velocity, for each spatial pixel a 
shock event is then defined to occur at the time when the velocity experiences an 
abrupt jump from large positive to large negative values (downward to upward 
motion). As a temporal threshold for this velocity jump, we utilize a maximum 
interval of two time steps, i.e. a 38 s window. The threshold for the amplitude of 
the velocity jump is set at a value of 4 km s . This is safely above the r.m.s value 
of ~1 km s -1 measured for chromospheric velocities at the dominant periodic- 
ity, and allows us to identify the most energetic events. Experiments with lower 
and higher values of the velocity threshold cause a different number of shock 
events to be counted, but do not alter the main conclusions about their overall 
spatio/temporal distribution. 

With these two combined criteria we can identify the presence of an acoustic 
shock at any given time or location within the FOV, and hence derive a temporal 
series of binary maps defining their spatial distribution, similar to the examples 
displayed in Fig. [5] a and e. A value of the shock strength can also be defined 
as the maximum velocity difference between downward and upward velocities 
during the shock. Given the limitations of the Doppler shift of the intensity min- 
imum, described above, both of these values can be considered as lower limits to 
the actual distributions. 



A.2. The Proper Orthogonal Decomposition approach 

As an independent validation of the results obtained through the velocity thresh- 
old method, we also attempted a more general approach, based on an analysis 
of the whole spectral shapes. This approach makes use of the Proper Orthogonal 
Decomposition (POD) technique, initially developed to study coherent structures 
in cha nnel flow turbulence and recently applied to various astrophysical prob- 
lems <Rees et aLteOOOklCarbone et all200alLawrence et alj2 004: Vecch io et alj 
2005; Vecchio 2006). As the POD technique is not yet widely known, we will 
describe it in some detail in the following. 



A.2.1. The POD method 

In the framework of the POD, a spatio-temporal field u(r, t) is decomposed as 
S %o a j(ty¥j( r )' The orthonormal basis functions, l F J , are not given a priori but 
obtained from the experiment, and can hence assume the proper functional shape 
of the phenomenon. The time coefficients at(t), computed from the projection of 
the data on the corresponding basis functions, represent the time evolution of 
the j-th mode associated with that eigenfunction. The modes are then ordered 
according to a parameter, the "energy" content of fluctuations associated with 
them, quantifying the relative contribution of each mode to the signal reconstruc- 
tion. The most energetical mode ( j = 0) is associated to the average contribution 
to the signal. For more details, we refer to the text bv lHolmes et al] 119961) . 



A.2.2. Application of POD to Ca II 854.2 profiles 

We apply the POD to analize the temporal evolution of the spectral profiles (like 
those shown in Fig. [4} for each pixel of the FOV separately, considering the 
(wavelength- dependent) intensity fluctuations around the mean profile. Thus, in 
our application the the POD decomposition of a wavelength-time spectrum will 
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be 

CO 

j=0 

namely, the role of spatial coordinate is played by the wavelength. For most of 
the pixels, we find that the 99% of the total energy is contained in the first four 
POD modes (including the average, j = 0), i.e. these are the only significant 
modes in the signal reconstruction. 

Panel e in Fig. |A.l| displays, as an example, the actual temporal evolution 
of the spectrum corresponding to an internetwork pixel in data set 1, where we 
find evidence for several shocks. Panel a shows the corresponding POD eigen- 
function *Pq(X) of the mode j = 0, i.e. the temporally averaged spectrum. The 
chromospheric steepening of acoustic oscillations, which eventually will lead to 
shock development, induce an intensity fluctuation characterized by a well de- 
fined behavior, i.e. a large displacement of the line from the average position, 
with an asymmetric contribution in the two wings. The POD analysis isolates 
this behavior in a single mode, whose eigenfunction has the typical shape shown 
in panel b of Fig. [TO] It is characterized by two opposite lobes of different ampli- 
tude, indicating positive and negative contribution to the average profile towards 
the blue and the red wavelengths. The sign of the time coefficient (to be multi- 
plied to the eigenfunction) selects when the positive contribution is in the red or 
in the blue part of the spectrum. 

Thus for each pixel the mode associated to the (possible) shock contribution 
can be identified by looking at the shape of the eigenfunctions, including the 
lobe asymmetry. Suppose that for the spectrum of Fig. IA. W e this is the mode 
j = i. The wavelength-time intensity contribution of this mode to the actual 
spectral evolution can be expressed as a,(?) l Pi(/i), and is shown in panel /of 
Fig. lA.lK for clarity we have displayed only its positive part). It illustrates how 
the modifications to the spectral profiles induced by this mode are in general 
stronger and more extended in the blue wing of the line, as well as the fact that, 
on average, the variations in the red wing last longer than in the blue wing. The 
latter property is present also in the simulations of CS97. For this particular pixel, 
£ = 1, i.e. the mode characterizing the eigenfunction ¥,(/t) is the one with the 
highest "energy" content of fluctuations. 



We further note that some spatial locations within the magnetic network 
areas display "shock- associated" eigenfunction £ = 1 as well, but are character- 
ized by longer periodicities (frequency below 4 mHz). They will be considered 
in more detail in a future paper. 



A.2.4. POD shocks 

Finally, in order to identify the occurrence of shocks from a pattern such as that 
displayed in Fig. lA.llf or h, we can simply count the events characterized by a red 
peak followed by a blue peak, with the jump occurring within the same temporal 
delay defined before (38 s), without applying any threshold to the amplitude of 
the events. Rather, an amplitude can be defined as the intensity difference of the 
two peaks of the eigenfunction (a, (r)T, (/i)) on either (temporal) side of the shock. 
Such difference is related to the variations induced by the shock over an average 
profile for the pixel, and is then again a measure of the strength of the event. As 
in the velocity analysis, repeating this calculation for each pixel we can derive 
a temporal series of binary maps, displaying the shocks' spatial distribution, as 
well as maps of their amplitudes. Panels a and e of Fig.|6]show typical "shocks 
maps", for a time around the middle of the temporal sequence. 



A.2.3. Dominant POD modes and periodicities 

Obviously, not all the pixels within the FOV have spectra similar to those dis- 
played in panel e of Fig. 1701 Panel g of the same Figure provides an illuminating 
example, showing the temporal evolution of the spectrum for a pixel pertaining 
to a fibril within the FOV of data set 1 . The average spectrum (POD eigenfunc- 
tion of the mode j = 0) is given in panel c, the "shock-mode" eigenfunction 
for j = i is shown in panel d, and the POD reconstruction for this mode is dis- 
played in panel h. For this case, the mode characterizing the eigenfunction ¥,(/!) 
is i = 3, i.e. with a much lower energy than the example of panel e. 

Panels a and b of Fig.fA/2]display the spatial distribution, for both data sets, 
of the mode £ defined above. (As said before, modes with £ > 3 contain an energy 
less than 1% and have not been considered.) For the majority of the internetwork 
pixels £ = 1, i.e. the "shock-mode" is the most energetic component of the in- 
tensity fluctuation. However, there exist large portions of the FOV, apparently 
pertaining to the photospheric internetwork but surrounding quite extensively 
the magnetic elements, for which the relevant eigenfunction corresponds to the 
mode i = 2. Their typical energy fluctuation content is less than half that pertain- 
ing to the pixels with shock-mode £ = 1, This effect is more prominent in data 
set 2, in correspondence to a larger extension of the fibrillar structures visible in 
the line core images. 

The coefficients adt) provide the temporal behavior of the different POD 
modes for each spatial pixel. For example, the coefficients used in the recostruc- 
tions of Fig. 1701 /' and /; are plotted in the same figure (panels £ and j), and show 
that for the first case the "shock" eigenfunction has a clear periodicity near- 
ing 150 s, while for the second one the temporal coefficient has a more erratic 
behavior. Maps of the dominant periodicities in the aj(t) coefficients, for each 
pixel and for the shock-associated j = i POD mode, are reported in panels c 
and d of Fig. |A.2| They have been derived from a Fourier analysis of the coef- 
ficients, but the same result is obtained if using a wavelet analysis of the time 
series, by considering the wave period containing the highest number of wave 
packets with significant power. In both data sets, for most of the pixels charac- 
terized by the £ = 1 POD mode, the dominant periodicity is ~ 200 s and slightly 
shorter (frequencies between 5 and 7 mHz). This is consistent with the typical 
chromospheric periodicity of ~ 3 minutes for acoustic oscillations, deriving from 
the atmospheric filtering of the photospheric p-modes. The association between 
i = 1 POD mode and this periodicity essentially reflects the strong steepening of 
the acoustic waves propagating in the chromosphere. We note that the two data 
sets display differences in this association, with a more fragmented distribution 
of the short-periodicity areas, and a stronger presence of longer periodicities (4 
mHz or less) in data set 2. 



18 



A. Vecchio et al.: The solar chromosphere at high resolution with IBIS 




10 III M H H t 11 tt » 40 » 

0**t d*ln> Mate f*itnl 

Fig. A.l. Example of POD analysis for two pixels from data set 1. Panels a-b, e-f refer to a pixel displaying a large number of 
shocks with 3 min periodicity, a: Eigenfunction of the mode j = (average of the signal), b: Eigenfunction describing the shocks 
contribution, displaying the typical asymmetric shape. Since the eigenfunction are orthonormal their amplitude is arbitrary, e: Actual 
measured spectra./: Contribution to the profiles in e reconstructed using only the POD mode associated to an eigenfunction like b 
(mode 1 in this case). For clarity, only the positive portion of the reconstructed function is displayed. Panels c-d, g-h: As before, 
but for a pixel where only few, energetically unimportant shocks are found by the POD analysis (shock mode = 3). Panels i—j: POD 
temporal coefficients used in the reconstructions of panels /and h. A3 minute periodicity is evident in panel i. 
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Fig. A.2. Panels a and b: maps of mode pertaining to the shock eigenfunctions, for data set 1 and 2. c and d: maps of dominant 
periodicity in coefficients a,(f), for data set 1 and 2. Spatial axes in arcsec. 



